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Abstract 



The boundary problem for a two-dimensional elliptical equation -sinh-Gordon 
has been investigated. The exact solutions have been found and identities of 
traces have been proposed. The application of the problem to the model of 
the Coulomb's plasma with the negative temperature has been considered. 



1. Introduction 



This paper is devoted to research and the construction of the explicit solutions of the 
nonlinear elliptical equation -sinh-Gordon: 



where A is two-dimensional Laplas's operator. 

On the one hand, it is of a great mathematical interest due to is describes the immer- 
sion of a negative curvature surface in the three-dimensional Euclidean space [1], and on 
the other hand it has a number of physical applications. In particular, it arises as the 
model of two-dimensional Colomb's system at the negative temperatures [2,3]. 

L.Onzager appears to have been first who introduced a conception of negative temper- 
atures in the context of a problem of a vortical line description. Further the conception 
turned out to be connected with so-called anomalous systems consisting of weakly interact- 
ing particles, spins and so on (see, for example [4]) where an energetic spectrum is bounded 
above and some maximal value of energy is available. In this case a statistical sum of 
the system is ^2 i g(E i )exp(—E i /T) under assumption that g{Ej) ~ exp(— aEi), E^ — > 
oo, a > 0, where g{Ej) is a density of a statements number , and Ei is an energetic level 
that can be finite both at T > and T < (T is absolute temperature). Assuming that 
a = oo and hence g(Ei) = at Ei > E m8/X , it is not difficult to see that the statistical sum 
changes into a finite series and has a finite value at T G (-co, oo). Then it is obvious that 
the density of the energetic levels that is proportional to Boltsmann's factor exp(—Ei/T), 
at T > drops as the energy increases. 

In particular, at T — > oo exp(— Ei/T) — > 1, i.e. all levels proved to be distributed 
uniformly. Another picture is at T G (— oo,0), in this case Boltsmann's factor is an 
increasing function of Ei and there exists an phenomenon of an inversion of a levels 
distribution with the density at the up levels becoming more than at down ones. It 
is clear, that the equation (1.1) involving the difference of Boltsmann's factors for the 
one-charged particles in its right side exactly corresponds to the situation of a negative 
temperatures that was described above. 

We now turn to a statement the boundary problem. Let us suppose that the function 



An = —4 sinh u, 




u(x,y) determined on the half-plane M. 2 ^ = {(x,y) : x G (—00,00), y > 0}, is the real- 
valued and enough smooth, u(x, 0), u y (x, 0) are boundary values of the desired initial 
function and its normal derivative correspondent ly. Also we assume that one of these 
condition is known, but there is some relation of the condition to one another that will 
by some given below (a nonlinear analog of the condition of the third kind). 
Moreover, we will assume, that 

u(x, 0) — > at \x\ — > 00. (1.2) 

Equation (1.1) is integrable and can be represented in the form of the compatibility 
condition of the linear matrix system: 

V x (x,y,\) = U(x,y,\)V(x,y,\), (1.3a) 

* y (x, y, A) = V(x, y, A)*(x, y, A). (1.36) 

Here ^{x,y, A) is a matrix-valued (2 x 2) function, A G C is a spectral parameter, U and 
V will look like: 

17(A) = U(x,y, A) = ^ + 2( A A ^ 1 1) (cosh M - l)]a 3 + 2 ( X + _\) sinh ue- 2ixa:i a 2 + 

V(X) = V(x,y, A) = +i ^T)( coshM " Qfa + |l3iy smh M e- 2 — a 2 + 

2 

It should be notice that a rather nontrivial and complex form of the gauge of the oper- 
ators U and V is connected with the gauge equivalence of the equation (1.1) and two- 
dimensional isotropic Heisenberg's ferromagnet, which was proved earlier [5,6]. The choice 
±1 as the poles is necessary to obtain a convenient system of contours in solving Riemann's 
problem (see below). 

2. Associated linear problem 

Let us ty ± (x, A) are the solutions of the (1.3) with the conditions (k(\) = 2/(A 2 — 1)): 

^ ± (x, A) -> exp(ik(X)xa 3 )(l +o(l)), x^±oo. (2.1) 
Let us determine a transition matrix T(A) = T(X,y) (A 2 = A 2 ): 

9-(x, A) = ¥+(*, A)T(A), T(A) = ( ^ ^ ) . (2.2) 

The matrix T(X) is unimodular: detT(A) = 1. 
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Assuming ± (x, A) = \l/ ± (x, A) exp(— ik(X)xa 3 ), from (1.3a) we obtain: 

0±(A) = e -" CT3 e (u/2),T1 ( T 2 e-" ,J3 ± (-A)cT 2 , A 2 = A 2 , (2.3) 

0±(A) = e - 2 ^a 1 ± (i)a 1) A 2 = A 2 . (2.4) 
From these relations it follows the relations of symmetry 

T(A) = <7 2 T(-A> 2 , A 2 = A 2 , (2.5) 

T(A) = o{T (1)^. (2.6) 

Then we have 

T(A) = (7 3 f(-i)(73, A 2 = A 2 . (2.7) 

We find a dependence of scattering data from the variable y. For this purpose let us 
turn to the system (1.3b). On fulfilling standard operations and taking /(A) = 2A/(A 2 — 
1) = Xk(X), we obtain: 

a(A, 0) = a(A, y), b(X, y) = b(X, 0)e 2 ^, (2.8) 

c(A, y) = c(A, 0)e- 2/ ^, d(A, 0) = d(X, y). (2.9) 

It follows from these relations that the coefficients a(A) and d(X) act as the produced 
functional of "the integrals of movement". Alternatively, from the requirement of finite- 
ness of the coefficients 6(A) and c(A) it follows, that it is necessary to put b(X, 0) equal to 
zero at 

A>1, -1<A<0, (2.10) 

and put c(A, 0) equal to zero at 

0<A<1, A<-1. (2.11) 

These inequalities determine boundaries of the zones. Domains in the real axis that are 
additional to these zones are continuous spectrum. There are no such restrictions at 
A = —A, so the continuous spectrum can be find in the imaginary axis. 
Taking U (x, y, A) = ik(X)a 3 + Q(x, y, A) , where 

/ |g(coshu-l) (-|f±isinhM+ f )e~ 2lx \ 
Q{x,y,X) - ^ (iA±i sinhu + ^ e 2ix _|A±i( coshM _ !) y 

and introducing the columns 4>f, 4>f of the matrix ± , from (1.3a) we will obtain Volterra's 
integrable equations of the direct problem: 

fc(x, A) = ei + J d^(x - X)Q(Z, A)0r(C A), (2.12) 
2 - (x, A) = e 2 + J d£g 2 (x - £, \)Q(£, A). (2.13) 
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In these equations t\ = (1, 0) T , e 2 = (0, 1) T and g 12 ( x , A) are Green's "initial" functions: 

g -{x, A) = 9(x)diag(l, e~ 2ikWx ), g^{x, A) = 9(x)diag(e 2tkWx , 1). (2.14) 
From (2.12)-(2.13) it follows that 

4>i(x, A) G H(£li U Q 3 ), 0j(a;,A) G ff(Q 2 Ufi 4 ), (2.15) 

where H(£li) is the class of functions (Qj = {A : (i — l)7r/2 < argA < iir/2},i = 1, ...4), 
that are analytical in the domain Qj. 

Completely analogously we can show that 

^(s,A)eff(n 2 ufi 4 ), $-(x,A) e if(fiiUfi 3 ). (2.16) 

From (2.12)-(2.13) we can easy obtain the integrable representations for the scattering 
data: 

«(A) = 1 + J <Wn(£, A)0 n (e, A) + Q 12 (£, A)0 21 (£, A)], (2.17) 
6(A) = | dflQnfo A)0r 2 (e, A) + Quit, A)0 22 (£, A)] e - 2 ^, (2.18) 

and hence 

a(A) G H(Qi U fi 3 ). (2.19) 

From (2.18) and the similar expression for c(A) it follows that generally speaking, b(\) 
and c(A) don't admit an analytical extension from both the real axis and imaginary one. 
On assuming y = 0, in (2.18), accordingly to (2.10), we obtain the equality: 

J ^[Qn(£, 0, A)0r 2 (£, 0, A) + Q 12 (£, 0, A)0 22 (£, 0, A)]e' 2 ^ = 0, (2.20) 

that connects the boundary values of the desired function and its normal derivative and 
acts as the boundary condition (in term of the spectral representation). 

Let us introduce a matrix complex- valued function Q = Q(x,y), that can be defined 

as 

Q(x,y) = lim V~(x,y,\). 

A|^oo 

From (2.21 ), (2.3) it follows, that the expression 

is fulfilled. This relation will be used later. 
Matrix Q obeys the equation: 

Q x = U(x,y, oo)Q, 

and, in addition, detQ(x,y) = 1. We suppose also that it has asymptotics 

Q = lim n(x) = (-l) N e^ a3 , 
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(2.21) 
(2.22) 

(2.23) 
(2.24) 



lim n(x) = I, (2.25) 

x — > — oo 

where cto is some parameter, a ^ 0, and the sense of the integer number N will be 
clear further. The choice of the asymptotics is in agreement with (1.3a), (2.23) and is 
stimulated by the gauge equivalence of boundary problems. 
From (2.2), (2.23), (2.24) we have: 

T(oo) = (-l) w exp(^ 3 ), T(0) = (-l) Jv exp(-^ O - 3 ), (2.26) 
and from (2.10) and the similar expression for c(A) it follows that 

6(±1 -0) = c(±l + 0) = 0. (2.27) 

3. Dispersion relations and identities of traces. 

The results are obtained in the previous sections allows to write down the dispersion 
relations and to obtain identities of traces as well. 
From (2.17) we have: 

a(X) = (-l) Ar e^ a [l + 0(^-)], |A| -> oo, AefiiUfl 3 . (3.1) 

A 

Similarly it can be shown 

d(\) = (-l) w e -^[l + 0(^)], |A|->oo, AGfi 2 Ufl 4 . (3.2) 

I I 

We assume now that the coefficient a(A) has finite number of simple zeros in domains 
Qi and f^, with the quantity of the zeros being same and equal to N in both domains. 
Let us also assume that the zeros don't belong to the continuous spectrum (may be 
except of the zones accordingly to the formula (2.10)). Moreover, we will consider that 
a(X) = a s (X)a c (X), where a s (X) is given by the contribution of the discrete part of the 
spectrum ("solitonic" part) and a c (X) - to the continuous part spectrum. Taking into 
account (2.2), (2.10)-(2.11) and equality |a s (A)| 2 = 1 at A = A, we will have 

lna c (A) + In rf c (A) = 0, A = A, 

lna c (A) + \nd c (X) = ln(l - |6(A)| 2 ), A = -A. 

These equalities assign the bound of piecewise analytical functions lna c (A) and lnd c (A) 
in the continuous spectrum. On applying the Cochy's formula to the function lna c (A) G 
H(Qi U Q 3 ) and taking into account both (2.26) and two previous relations, we obtain 

, f°° dfi ln(l - \b(ifi)\ 2 ) . . , 

In a c (X) = - / -^-5= ' V WI > sign(fi). (3.3) 

2m jj, + i\ 



CO 



Now let us consider the "soliton" part of the coefficients a(A): a s (A). Let a s (X n ) = 
0, where A n G Q±. Then in a view of the involution a(X n ) = a(— 1/X n )e ia ° we have: 
a(— 1/A) = 0, — A -1 G ^3. From this and with consideration for (2.6) we will find the 
expression for a s (X): 
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) ' 



(3.4) 



Combining the equation (3.3) and (3.4), we will finally obtain 





(3.5) 



From the latter with consideration for (2.26) it follows: 




(3.6) 



e 



— ICtO 




(3.7) 



It follows here from that at a = and N — 1 either Ai = Ai or Ai = — Ai. In the 
first case the eigenfunction proves to be non-localized and the solution doesn't fit to the 
asymptotic (2.1). In the latter case (Ai = — Ai) a s (Xi) = 1; means that discrete eigenvalue 
is in the continuous spectrum (the situation that is beyond the consideration). Hence in 
the reflect ionless case at = the minimal value of N is equal two. This explains the 
cause of introducing the parameter a in (2.24). 

Assuming A n = p„e ie ™, p n > 0, 9 n e (0, 7r/2), in (3.7), we will obtain Ylu=i^n — 
—a /4: + nn/2, n — 0, ±1, ±2, ... ; in particular, at N — 1: Q\ = —ao/4 + nk/2, where k 
is the number to fulfill the inequalities < — «o/4 + nk/2 < ir/2. 

Let a s (X) = 1, i.e. discrete spectrum in the system is unavailable. In this case from 
(3.6) we have: 



It should be noticed, that the expression of the type of (3.5) holds for the coefficient d(X) 
also, however we will not discuss it here. 

Let us turn to deducing the identities of traces. Do to the matrix Q (the relation 
(2.23)) can't be determined explicitly the equalities in the points zero and infinity can't 
be obtained (because we don't know solutions of the Riccati's equation in this cases). 
Therefore we will restrict ourselves by the consideration of the relations near points of 
the poles divisors. 




(3.8) 



From (3.5) and (3.8) it follow, that 



T(±l) = 1 . 



(3.9) 
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Assuming lna(A) = J2^L a^^iX ± l) p , A — > ±1, from (3.5) we have: = —dp *\ 

where p is a even number, and dp 1 " 1 = a p 1 \ where p is an odd number. These equalities 
being the consequence of the spectral structure of the problem permit to control the 
validity of intermediate calculations and the final result. 

From (1.3a) using the standard way of calculations one can obtain Riccati's equation 
for the function F(x, A) = 0^(x, A)/0j" 1 (x, A): 

F x + Q 12 F 2 = (Q 22 - Q n - 2ik(X))F + Q 21 . (3.10) 
Also from the equation for the function in view of (2.2), it follows that 

/oo 
dx[Q n (x, A) + Q 12 (x, X)F(x, A)]. (3.11) 
-oo 

Using explicit form of the matrix elements Qij, at the vicinity of the point A = — 1 we will 

obtain (F(x, A) = E^i^'^O + 1 Y))- F i~ 1] = (iu z /2)e 2ix . Then accordingly to (3.11) 
lna(A) = 0, A — > —1, which is in agreement with (3.9). In the next order (~ (A + 1)), 
considering (3.10), (3.11) and (3.5), we will have: 

Pj-l [°° dp (p 2 -l)\n(l-\b(ip)\ 2 ) 

+ p\ + 2 Pl cos X + 1 + 2tt _ i)2 + V " 

i f 00 

= — / dxu x (x, 0)u y (x, 0), (3.12) 

J -co 

pi sin Z" 00 d////ln(l - |6(i/x)| 2 ) 



signal) = 



p\ + 2 Pl cos 6 1 + 1 2tt (p 2 - l) 2 + 4/i 2 

i r°° i 

= --y ^[(cosh W (x,0)-l)--(^(x,0)-^(a:,0))], (3.13) 

where Ai = pie tdl , Q\ e (0, 7r/2). Thus, the continuous spectrum of the problem con- 
tributes the left parts of the identities of traces. It should be noted that the similar 
character the identities have also in the neighbourhood of the point A = 1. 



4. Equations of the inverse problem and formulas for reconstruction. 



It follows from the analysis of properties of the operator spectrum of an associated 
linear problem and analytical properties of Jost's solutions performed above that the 
continuous spectrum of the problem is in "the cross": {A : Ar = U Xj = 0}, whereas 
matrix columns of these solutions admit the analytical extension to certain domains of the 
plane of a complex spectral parameter A. This allows using Riemann's technique problem. 
However there exists some additional difficulty here: from (2.21) and the determination 
of the matrix _ it follows, that <j>~(x, oo) = ^(x, oo) = Q(x,y), with the matrix Q(x, y) 
being unknown. This leads to Riemann's local problem on an analytical factorization 
appears to be characterized by noncanonical conditions at the infinity in this case. Also 
some difficulties arise in searching for formulas of the reconstruction. This section is 
devoted to these problems. 



7 



As the method of solution of the inverse problem we choose the method of Riemann's 
vector problem. From (2.2) and the connection ^ f± and (j)^ we obtain: 



<j ) -(x,\)=<j ) + (x,\)T l (x,\), (4.1) 

where 

Ti(x, A) = exp(ik(X)xa 3 )T(\) exp(—ik(\)xa 3 ). 
We use one of the equalities following from (4.1) 

^ = «M)-^^. (4.2) 

Now let us suppose that A n G fix is such that a(X n ) = and also — 1/A n G and 
a(— 1/A n ) = 0. Then the vector-functions ^(x, A) and ^(x, A) proved to be linear- 
dependent at these points, from here we have: 

07(x, A n ) = --J-0+(x, A n )e- 2ifc ( A »)*, (4.3) 

ft(x, -f) = i^(x, - J-) e -**<-£>*, (4.4) 

where 0^(A„), 0j(A„), 0r(— A^ 1 ), 0^(— A^ 1 ) are eigenfunctions of the discrete spectrum, 
b n are the coefficients of transition. 

Taking into account that (A) — > f^(x, y), |A| — > 00, A G f2iUf23, and <fti (X) / d(\) — > 
f2}(x, y), A G f^Uf^, where f2} = Q(x,y)ei, and applying Cochy's formulas to the vector- 
column (<f>i—Q\) that is analytical in fli and also using (4.2)-(4.4), we obtain the following 
equality: 

N R 1 

0r(A) = n\ - iY, A2(; f% ) e- 2 ^ ) -x(x)0r(-f)+ 

~[A + A n y_oo 27T a(/i) /i + A 

- /" ^^e--<^rM« S »M, (4-5) 
J_ 00 2 7r a(z//J /J, + i\ 

where x(^) = e - ixa3 a 2 e~ iul2ai e- ixa \ (3 n (y) = b n /a'(X n ) = (3 n (0)e 2l ^ y . The equation 
is not closed equation for the column <fii(\) yet due to it involves an unknown vector- 
function fl\. Therefore to obtain the necessary equation we will use (2.22) rewritten in 
the form: 

n- l e- ix ^a 2 e-^ 1 e- xa m = a 2 . (4.6) 
Substituting (4.6) into (4.5) and putting 

0r(A)=Q(x,y)0r(A), (4.7) 
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we obtain the following system of singular integral equations for the function ^ (A): 

N - p -2ik{\ n )x 

A + A„ 



0i ( A ) = e i + i Yl ^ n6 \ T T ( A «) 



n=l 



Ll^ rT a 20i(-^)+ / x - T i °-20i {n)sign(n)- 

^t A n( A -A^) A ™ 2?r a(/i) /i + A 



Z" 00 d/x 6(i/x) e ~ 2ik ^ x 



N 



£j^ e -2ik(\n)x 



0i ( A p) = ei + i V " = — ^20i (A n )- 

n=l A P + A « 

. A p n e~ 2ik{ ^ )x ~_. 1, f°° dfib(fi) e- 2lk ^ x ~_, . . 

n^ A n( A p-A^) A « J-oo 2?T a(/j) /J + A p 

- f _ u .. aiKMsignM, (4.9) 
1 f3 n e~ 2ik ^ x 

0r(-^) = e i + i z2 t — 77^ (j 20r( A n)+ 

A p n=1 A n - 1/ A p 

, .^ A^ 2 *^ -- 1 d/x 6(/z) e- 2 - fc ^- = 

t^i X "(x- P + XT) A « ^-oo 2tt a(/i) /i - ^ 

dfibiifi) e~ 2ik ^ x ~- , , , , 

-^0-201 {i[J,)sign(iJ). (4.10) 



The equations (4.9), (4.10) are obtained from (4.8) at A = A p and A = — 1/A P respec- 
tively. On the right-hand sides of these equations it is necessary to take into account the 
availability of the zones (2.10). 

The system (4.8)-(4.10) is the system of the equations for the eigenfunctions of the 
continuous and discrete spectrum. Considering the system is written for the "gauge" 
vector-functions, the necessary formulas of the reconstruction should be written in terms 
of the same function. 

Using (1.3a), (4.7), we obtain the differential equation for the matrix 0~ = <f>~(x, A): 

0- = -ik(\)4>-<j 3 + — !— ft _1 [i(coshu - 1)<7 3 + sinhwe- 2ixa3 c72]fi0~ + 
A — 1 

+ik(\)n- 1 a 3 ncf)-. (4.ii) 

In the neighborhoods of the point A = 1, setting 0~ = X]fc=o°0 fc ( A ~~ (0o — * x ~ ¥ 
— oo), we find: 

^-'(tanhW" t<mh f" 2 "). (4.12) 
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Now we parameterize the unknown matrix f2 by four complex functions of variables 

x, y: 



Q(x,y)=[ a ^ y \ % X >V\ ). (4.13) 



l(x,y) 5(x,y) 
Then the relation (4.6) gives 4 algebraic connections: 

a = 5 cosh — — 3 sinh — e~ 2iX , 3 = a sinh — e ~ 2lx — 7 C osh — , 
2 2 2 2 

7 = -j3 sinh ^ + 5 cosh ^-e 2tx , 5 = a cosh ^ - 7 sinh ^e 2 " 1 . 

Z Z Z Z 

In a view of the relation and from (4.12) we obtain (0° = (f) ): 

^ - Jrf T - < 4 - 14 > 

where det(/>° = (cosh|)~ 2 . Combining (4.12), (4.14) and (4.6), we will have: 

cosh ^ = [|0? x | 2 + l^l 2 ]- 1 = [|0? 2 | 2 + l^l 2 ]" 1 , (4.15) 

or 

cosh^!^ = [< > 2 ]- 1 = [< 0~f ,0° > 2 ]-\ (4.16) 

where the symbol <, >2 denotes the scalar product of the vectors in C 2 . Another equiv- 
alent representation of (4.16) is 

cosh^%^ = {l-Tr^} )}- 1 . (4.17) 

Zj z 

Analogously we can obtain one more formulas of the reconstruction: 

, u(x,y) 2Ae~ 2ix 
sinh = , (4.18) 

2 1 + V 1 - 4A 2 e~ 4ix 

where A =< (((0~)°)i) T , (((</>~))°)2 >2, A — > at \x\ — > 00, and the requirement of reality 
of the potential leads to the value Ae~ 2lx is to be real. It is easy to prove the latter using 

(4.14) and the definition of A. Furthermore, employing elementary transformations it is 
demonstrated that (4.18) is equivalent to (4.16). 

From (4.15), (4.18) it follows that the solution of the problem (1.1), (1.2) has a non- 
singular character. 

In solitonic sector of the model (6(A) = 0) one can be found the expression for the 
formulas of the reconstruction through eigenfunction of the discreet spectrum, can find 
the expression in terms of the eigenfunctions of discrete spectrum. For this purpose we 
first resolve the right of (4.8) at A — > 1 and after equating it to $? we will obtain from 

(4.15) the following: 





2 


+ 


rj 


2^2 


iei 2 + 


7/ 


2 



u(x, y) = 2 In v „ nJ , (4.19) 
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where 

2N j3 n e~ 2ik( ^ n)x ~ 

£ = y) = J2 n \ ii — ^n( A «)' 
n=i K + 1 



2N 



e -2ik(\„)x - 



v = v (x,y) = l + ^2 n - — 4> 2i (A„), 



An+JV — — T^, Ai+JV — — "p^, 71 — 1, 2, ...A?". 
•*n A r 



The functions n (x, A„), 21 (x, A„), involved in (4.19) can be found from the system of 
linear algebraic equations (4.9)-(4.11). Then: 

f =< AE, R^E > 2n , rj = 1- < RE, R~ l AE > 2n , (4.20) 

where: 

^4 — {A mri }, A mn = E m E n / {\ m -\- A n ), -£/m = /3m 6 ^ m ^ , 

A = ^a( ? ((A 1 + l)- 1 ,....(A 2n + l)- 1 ), m, n = l,....2N, E = (E u ....E 2N ), 

R = I + AA. 

N-"soliton" solution (4.19), (4.20) is the exact solution of the initial boundary problem. 
However in the simplest case (N = 1) it is it is rather cumbersome and it is beyond the 
paper. 

Also it should be noted that the asymptotics of the solution for the continuous spec- 
trum has been calculated in [6] proceeding from gauge equivalence between the given 
problem and the problem for Heisenberg's ferromagnet, with the technique of the ref. [7] 
being used essentially. 



5. Application to the model of one-charged plasma 
at negative temperatures. 



This model is lead by the model of the equilibrium system of charged particle (vortices) 
on the plane [2, 3]. The latter permits generalization of Debay-Hukkel classical model [8] 
for plasma or an electrolyte solution to the case of a nontrivial statistical "phone". 

Following [2], [3], let us suppose that there is the Coulomb's system of N$ particles 
possessing a charge +e and Nq particles possessing a charge— e embedded in some finite 
volume V = L 2 (where L is the linear scale). Furthermore, we will suppose that a 
temperatures of the ionic component coincides with a temperature of electronic one, with 
both of them being same and equal to T (Boltsmann's constant is considered to be equal to 
unit). Using the BBGKI chain it is convenient to consider a scalar potential $ = y): 

e$(x,y) = n J dx'dy'<f)(r - r')p(r'), r=(z,y), (5.1) 

where Uq is the average density of charge particles, 0(r — r') is Coulomb's two-particles 
interaction, p(r) is a difference between the ionic density and an electron one. Then 
employing rather rigorous procedure of the reduction of the chain we can obtain the 
equation: 
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A$ = 



[exp(— ) - exp( 



~T 



)] 



(5.2) 



where 1° is a characteristic length, 



V 



J dxdy[exp[^ J dx'dy'(j)(r - r')p(r')]] 
Setting u = <f and Debay's radius A = (AX 2 D /c)A, X 2 D = Z°|T|/(87rn e 2 ) in (5.2), we find: 



This equation leads to the equation (1.1) at T < ( the case T > is considered in 



According to the statement of the problem describe above the equation (5.3) held for 

should be added by the given value of the field at the boundary considering. 

Now we consider some physical aspects of the results obtained in the previous sections. 
Assuming n = Nq /V = const, i.e in the thermodynamical limit. 

For this purpose we will reveal some features of the behaviour of the equation (4.19). 
The total energy of the finite system of discrete isolated charges for the two-dimensional 
case [2], [3] is: 



where iV + = £\ iV +, JV " = £\ Nq., $ - = -(2e 2 // ) In |n — r j | is Coulomb's potential of 
the interaction between the particles of "i" and "j" clusters, r ; and rj are their radius- 
vectors. 

Comparing the expressions for $jj with (4.19), we can see they have the similar form. 
Let B = |£| 2 + |?7| 2 , and < B < 1, due to the potential is real. Then from (4.19) it 
follows, that u(x, y) ~ 2 ln(f?/2) — > — oo at B — > and u(x, y) — > at B — > 1. Here in the 
terms of eigenfunctions of the discrete spectrum and values £, 77, the first condition means 
that: ^, 77 — > 0, which leads to ^ 2 l ^ 1 [ / 9„exp(— 2i/c(A n )x)/(A„ + l)]0^ 1 (A n ) — > —1; it is 
possible on the case of finite r. The latter can be fulfilled, in particular, at £ — > 0, 77 — > 1, 
i.e. at r — > 00. 

Thus, we have been obtained the following qualitative picture. Let us suppose there 
is no external field first {u(x, 0) = u y {x, 0) = 0). Also clusters of the same charges can 
by formed from the charged particle. After appearance of an external field the space 
distribution being available before changes into another stationary state. At the any 
point (x, y) G M 2 ^ the potential is the sum of potential arisen by the the boundary dM. 2 ^ 
and potential of system; as result, some self-consistent field is created. The potential of 
the field is described by expression Tu(x,y)/e, where u(x,y) = u(B), and the function 
B = B(x, y) plays the role of the "transmission" function of the medium. At certain finite 
x, y a local and deep minimum of the potential corresponds to the minimal value of B, 
with u < at the vicinity of the minimum. Considering T < a value $ proves to to be 
big and positive. This fact shows a local concentration of particles of the identical sign. 
Due to non-singular character of (4.19) $ keep be bounded up and it leads to a limitation 
of the energetic spectrum of the system. So we deal with an abnormal system where 



sign(T)Au = 4sinhw. 



(5.3) 



[9, 10]). 




(5.4) 
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the inversion of the level distribution is available, it means that the positive value of the 
energy together with negative temperatures leads to the higher levels (in average) become 
more preferred for the system than the lower ones. Furthermore, the charges distribution 
of clusters with same sign appears more apparently. 

The condition r — > oo and u(B) — > corresponds to the maximum of B, i.e. as we 
recede from the boundary the field created by those charges fads due to the interaction 
with own field of the medium. The damping is slower than in the similar situation when 
temperatures are positive [9,10]. This is not difficult to see on comparing the solution of 
corresponding linear problems: 

Au = ±4w, u — > at r — > oo, (5.5), 

where the sign "+" corresponds to system with T > 0, whereas "-" corresponds to negative 
temperatures. In the former case Green's function is proportional to McDonald's function 
of the zero index and has the asymptotic: ~ (e~ x / \/x)[l + 0(1/ x)], x — > oo, and in 
the latter case it is proportional to Neemann's function of zero index that behaves like: 
(sin (a; — 7r/4)/y / a : )[l + 0(1/ x)], x — > oo. Therefore the slower damping of the field at 
T < seems to be explained by weaker defense of the boundary due to an alternation of 
domains with the same charged particles. 

The picture describe above is true everywhere on the semi-plane except of the certain 
directions defined evidently. In these cases B takes an intermediate values: Be [e, 1 — 
e], e > 0, and the potential is an oscillating function in those directions, i.e. there is the 
distant order in the medium. It means that a coherent structure of the distribution of an 
electric field and densities of a plasma distribution arises. 

The author is grateful to V.D.Lipovski for useful discussions. 
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